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Abstract. Heavy-ion reactions especially those induced by radioactive beams provide useful information about the density 
dependence of the nuclear symmetry energy, thus the Equation of State of neutron-rich nuclear matter, relevant for many 
astrophysical studies. The latest developments in constraining the symmetry energy at both sub- and supra-saturation densities 
from analyses of the isopsin diffusion and the 7t^ /k^ ratio in heavy-ion collisions using the IBUU04 transport model are 
discussed. Astrophysical ramifications of the partially constrained symmetry energy on properties of neutron star crusts, 
gravitational waves emitted by deformed pulsars and the w-mode oscillations of neutron stars are presented briefly. 
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CONSTRAINING THE DENSITY DEPENDENCE OF NUCLEAR SYMMETRY 
ENERGY WITH HEAVY-ION COLLISIONS 

To determine the density dependence of the nuclear symmetry energy Esym{p), thus the equation of state (EOS) of 
neutron-rich nuclear matter, has been a longstanding goal of both nuclear physics and astrophysics. The Esym{p) is 
critical for understanding not only the structure of rare isotopes and heavy-ion reactions IHHSHIlll, but also many 
interesting issues in astrophysics |@, 0, HI 01 ■ In this contribution we first summarize some recent progress in con- 
straining the Esym{p) from analyzing the isospin diffusion and 1%'^ ratio in heavy-ion collisions within an isospin 
and momentum dependent transport model IBUU04 IJ0i1 . We will then discuss some astrophysical ramifications of 
the partially constrained Egymip)- Some interesting information about the Esym{p) has been obtained over the last few 
years from heavy-ion collisions I's']. As an example, shown in the middle of Fig.[T]are the (p) used in the IBUU04 
transport model. The three MDI Esym{p) 111] are obtained using a modified Gogny force by adjusting the parameters 
introduced in the interaction. For comparisons, the IQMD 1IT2I1 and the APR llTsIl predictions are also shown. Shown in 

the left window is the degree of isospin diffusion 1 — Ri versus the slope parameter L = 3po — ^sp"^"^ \p=Po 

of the MDI symmetry energy. Within the IBUU04 model analysis, the MSU data of Tsang et al. fT?] favors a (p) 
between that with x = and x = —I around the APR prediction. More quantitatively, one can infer that the slope 
parameter is about L = 88 ± 25 MeV 1 17]. We notice that a very recent re-analysis of the MSU isospin diffusion data 
using the ImQMD model found that the extracted range of the (p) overlaps with that from the IBUU04 analysis 
around the E^mip) with x = (01 APR) 1 18]. In the MSU isospin diffusion reactions the maximum density reached is 
about 1 .2po ilSn. Moreover, the isospin diffusion was found to probe the symmetry energy during the expansion phase 
of the reactiori lM \l5t] . Thus the isospin diffusion data provides us information about the Esym (p ) in the sub-saturation 
density region. 

Very interestingly, a recent IBUU04 transport model analysis of the FOPI/GSI /n^ ratio data 1^2^ in relativistic 
heavy-ion collisions at SIS/GSI with beam energies above 400 MeV/A provides circumstantial evidence suggesting 
a rather soft nuclear symmetry energy at p > 2po compared to the APR prediction 1 19]. Shown in the right window 
are the calculated %^ 1%^ ratios in comparison with the FOPI data at 0.4 AGeV with the reduced impact parameter 




FIGURE 1. Left: Degree of the isospin diffusion 1 — Ri as a function of the slope parameter L of the symmetry energy. Middle: 
The density dependence of the symmetry energy. Right: The /k^ ratio as a function of the neutron/proton ratio of the reaction 
system at 0.4 AGeV with the reduced impact parameter of b/bi„ax < 0.15. The inset is the impact parameter dependence of the 
7t'/n+ ratio for the ^^Ru+'^^'Ru reaction at 0.4 AGeV. Taken from Refs. fTvlfTgll. 

b() = b/b„,ax < 0.15 as a function of the neutron/proton ratio of the reaction system. The inset shows the /n^ ratio 
as a function of bo for the ^^Ru+^^Ru reaction at 0.4 AGeV. It is seen that both the data and the calculations exhibit 
very weak b^ dependence for the / ratio, even for mid-central reactions where we found that the multiplicities of 
both and vary appreciably with the b^. For the symmetric '*''Ca+*''Ca and the slightly asymmetric ^^Ru+'^^'Ru 
reactions, calculations using both x — I and x = can well reproduce the FOPI data. Most interestingly, for the more 
neutron-rich reactions of ^^Zr-n^^Zr and '^^AuH-'^^Au calculations with x = 0,0.5 and 1 are clearly separated from 
each other. The FOPI data favors clearly the calculation with x = 1 . The corresponding symmetry energy is much 
softer than the APR prediction. A detailed study on the excitation function of the %^ jn'^ ratio also indicates that 
the FOPI data can be approximately reproduced only with the x=\ symmetry energy [19]. It is also interesting to 
mention that IQMD calculations 112011 give similar results as the IBUU04 with x = 0. This is not surprising since the 
symmetry energy functionals used in the IQMD and the IBUU04 with x=0 are very similar for po < p < 3po as shown 
in the middle of Fig. [T] The most importance influence of a rather soft symmetry energy at supra-saturation densities, 
such as that with x= 1, on pion production is through the rather high neutron/proton ratio reached in the participant 
region due to the dynamical isospin fractionation. As shown in the middle window, with x~ \ the Esy,„{p) at p > 2po 
reached in the reaction is very small. Thus a rather high N/Z in the participant region is energetically favored due to 
the dynamical isospin fractionation ll2lll22ll and thus the larger /n^ ratio is observed. In the reactions considered 
at 400 MeV/A, the maximum central density reached is about 2.5po. By varying separately the Esyi„{p) at sub- and 
supra-saturation densities used in the IBUU04 simulations for these reactions it was found that the /n^ ratio is 
much more sensitive to the variation of the Esyi„{p) at supra-saturation rather than sub-saturation densities. 

Putting together the information from analyzing both the isospin diffusion and the / ratio data, we expect that 
the Esyinip) reaches a maximum somewhere between po and 2po before it starts decreasing at higher densities. This 
indicates the importance of mapping out the Esym{p) at densities from about po to 2.5po. Such experiments are being 
planned at several facilities. If the Esym{p) is confirmed by more experimental and theoretical studies to decrease with 
increasing density above certain density, it not only posts a serious challenge to some nuclear many-body theories but 
also has important implications on several critical issues in nuclear astrophysics, such as, the cooling of proto-neutron 
stars, the possible formation of polarons due to the isospin separation instability |23, 24], the possible formation of 
quark droplets |25] and hyperons |26] in the core of neutron stars. 

Since the Esyi„{p) is only partially constrained in some density regions by the available experimental data within the 
IBUU04 transport model analyses, astrophysical applications of these constraints involve some interpolations and/or 
extrapolations under some assumptions. We notice here, however, that the extrapolation of any low-density symmetry 
energy to supra-saturation densities can be very dangerous as illustrated in Fig. 1. Using the same IBUU04 transport 
model, at sub-saturation densities the isospin diffusion analyses favors clearly the APR prediction, but above about 
2po the %^ 1%'^ ratio data favors instead the Gogny prediction that is much softer than the APR prediction. Fortunately, 
while many astrophysical questions depend on the symmetry energy over the whole density range, some issues, such 
as the core-crust transition density in neutron stars depends only on the Esym{p) at sub-saturation densities. In the 
following section, we present several examples illustrating the astrophysical importance of the (p) using the MDI 



EOS with X = and x— — I extrapolated to supra- saturation densities. Effects of the softening of the E^ym (p ), such as 
that with X = 1, is currently under investigation. 



NUCLEAR CONSTRAINTS ON PROPERTIES OF NEUTRON STAR CRUSTS 

Neutron stars are expected to have a solid inner crust surrounding a liquid core. Knowled ge on pr operties of the crust 
plays an important role in understanding many astrophysical observations ll6ll7i Isl I9I IztI l28l I29I Isol [31I [32[ [331 [34ll . 
The inner crust spans the region from the neutron drip-out point to the inner edge separating the solid crust from 
the homogeneous liquid core. While the neutron drip-out density Po„, is relatively well determined to be about 
4 X 10'^ g/cm^ [35], the transition density p, at the inner edge is still largely uncertain mainly because of our very 
limited knowledge on the EOS, especially the density dependence of the symmetry energy, of neutron-rich nucleonic 
matter . Recently, using the equation of state for neutron-rich nuclear matter constrained by the isospin diffusion 
data from heavy-ion reactions in the same sub-saturation density range as the neutron star crust, Xu et al. put a tight 
constraint on the location of the inner edge of neutron star crusts Ii36l] . 



The density and pressure at the core-crust transition 



The inner edge corresponds to the phase transition from the homogeneous matter at high densities to the inhomo- 
geneous matter at low densities. In principle, the inner edge can be located by comparing in detail relevant properties 
of the nonuniform solid crust and the uniform liquid core mainly consisting of neutrons, protons and electrons {npe 
matter). However, this is practically very difficult since the inner crust may contain nuclei having very complicated 
geometries, usually known as the 'nuclear pasta' |7, 32, 37, 38, 39]. Furthermore, the core-crust transition is thought 
to be a very weak first-order phase transition and model calculations lead to very small density discontinuities at the 
transition Il30ll40ll4lll42ll . In practice, therefore, a good approximation is to search for the density at which the uni- 
form liquid first becomes unstable against small amplitude density fluctuations with clusterization. This approximation 
has been shown to produce very small error for the actual core-crust transition density and it would yield the exact 
transition densit y fo r a second-order phase transition ifsol l40l l4li l42| . Several such methods including the dynamical 
method J^ypllH li^ IHB, the thermodynamical mefliodfll, IH IH] and the random phase approximation 
(RPA) 1I42II47II have been applied extensively in the Uterature. 
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FIGURE 2. Left: The transition density p, as a function of L (left windows) and if,,,,, (right windows) by using both the dynamical 
and thermodynamical methods with the full EOS and its parabolic approximation. Right: The transition pressure P, as a function of 
Pt and L within the dynamical method with the full EOS and its parabolic approximation. The MDI (upper windows) and Skyrme 
interactions (lower windows) are used. Taken from Ref. |36]. 

Shown in the left panel of Fig.|2]is the pi as a function of L and Ksym = ^Po~~^^T^\p=po using both the dynamical 
and thermodynamical methods with the full EOS and its parabolic approximation (PA) from the MDI interaction with 
the varying x parameter and 47 Skyrme forces ll36ll . With the full MDI EOS, it is clearly seen that the Pt decreases 
almost linearly with increasing L for both methods. This feature is consistent with the RPA results and that found 
recently by Oyamatsu et al. ll43ll . The similar relation is also observed between the p, and Ksym due to the fact that 
Ksym always correlate with L for a fixed energy density functional i36ll . It is interesting to see that both the dynamical 




FIGURE 3. Pt as a function of p, by using the dynamical method without parabolic approximation for both MDI interaction and 
SHF calculations. The shaded band represent the constraint from the isospin diffusion data. Taken from Ref. 1361 . 



and thermodynamical methods give very similar results. On the other hand, surprisingly, the PA drastically changes 
the results, especially for stiffer symmetry energies (larger L values). Also included in the left panel of Fig. |2]are the 
predictions by Kubis using the PA of the MDI EOS in the thermodynamical approach |45]. The large error introduced 
by the PA is understandable since the )3 -stable npe matter is usually highly neutron-rich and the contribution from 
the higher order terms in 5 is appreciable. This is especially the case for the stiffer symmetry energy which generally 
leads to a more neutron-rich npe matter at subsaturation densities. In addition, simply because of the energy curvatures 
involved in the stability conditions, the contributions from higher order terms in the EOS are multiplied by a larger 
factor than the quadratic term. These features agree with the early finding that the pt is very sensitive to the fine 
details of the nuclear EOS. Applying the experimentally constrained L values of 86 ± 25 MeV to the Pt~L correlation 
obtained using the full EOS within the dynamical method shown in the left panel of Fig. |2l one can conclude that the 
transition density is between 0.040 fm^^ and 0.065 fm^^. This constrained range is significantly below the fiducial 
value of p, « 0.08 fm^^* often used in the literature and the estimate of 0.5 < pt/po < 0.7 made in ref. |8] within the 
thermodynamical approach using the parabolic approximation of the EOS. 

The pressure at the inner edge is an important quantity related directly with the crustal fraction of the moment of 
inertia which can be measurable indirectly from observations of pulsar glitches It is very instructive to quote the 
analytical estimation obtained by Lattimer and Prakash |8fl for the transition pressure 
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where Kq is the incompressibility of symmetric nuclear matter at po and 5, is the isospin asymmetry at p,. One can 
see that, besides the implicit dependence on the symmetry energy through the p, and 5,, the P, also depends explicitly 
on the value and slope of the Esym{p) at p,. Thus the P, depends very sensitively on the £'.„„,(p). Shown in the right 
panel of Fig. |2]is the Pt as a function of Pt (left windows) and L (right windows) by using the dynamical method with 
and without the parabolic approximation. The results from Eq. ([TJ using the p, and £„.,„ coiTesponding to the full EOS 
and its PA are also shown for comparisons. It is interesting to see that the Eq. ([T]i predicts qualitatively the same but 
quantitatively slightly higher values compared to the original expressions for the pressure with or without the PA for 
both the thermodynamical and dynamical methods even though this formula was derived from the thermodynamical 
method using the PA. This observation implies that the direct effect of using the full EOS or its PA on the pressure 
is small although the PA may affect strongly the transition pressure P, by changing the transition density p,. The Pt 
essentially increases with the increasing p, in calculations using the full EOS, but a complex relation between the P, 
and Pt is obtained using the PA. The large difference in P, is due to the strong PA effect on the p,. Moreover, the latter 
does not vary monotonically with L for the PA as shown in the right panel of Fig.|2l 

It is also interesting to examine the range of P, corresponding to the p, and L constrained by the heavy-ion reaction 
data. In Fig. [3] we show the P, as a function of p, by using the dynamical method and the full EOS for both the MDI 
(solid line) and the Skyrme (filled squares) calculations. It is interesting to see that the MDI and Skyrme interactions 
give generally quite consistent results. Corresponding to the p, constrained in between 0.040 fm^^ and 0.065 fm^^, 
the Pt is limited between 0.01 MeV/fm^ and 0.26 MeV/fm^ with the MDI interaction as indicated by the shaded area. 
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FIGURE 4. Left: The whole radius R, the crust thickness hR, the core radius R, as functions of L at fixed total mass of 0.4Mq, 
1 .OM,;, and 1 .4M0 , respectively. Right: The crustal fraction of neutron mass AM /M, the moment of inertia / of the whole star and 
the crust contribution A/ as a function of L, at fixed total mass 'dAMr.,, \ S)M,-., and \AMr.,, respectively. Taken from Ref. |36fl. 



which is significantly less than the fiducial value of P, k, 0.65 MeV/fm^ often used in the literature [SJ. As pointed 
out in a recent work by Avancini et al [49], the value of « 0.65 MeV/fm^ may be too large for most mean-field 
calculations without the PA. We note that among the 47 Skyrme interactions used here, the following 8 interactions, 
i.e., the SkMP, SKO, R<j, Ga, SV, SkI2, SkI3, and SkI5, are consistent with the constraints from heavy-ion reactions. 
These results indicate that one may introduce a huge error by assuming a priori that the EOS is parabolic for a given 
interaction in calculating the p, and P, . 



The core size and crust thickness 

For a given neutron star of total mass M and radius R, what are the respective sizes of its core and crust? How do 
they depend on the stiffness of the symmetry energy? How do they depend on the neutron star mass M? How does 
the thickness of neutron star crusts depend on L while it is well known that the size of neutron skin in heavy nuclei 
increases with the increasing L 1 17]? These questions have been investigated recently in Ref. |36]. 

Show in left panel of Fig. |4]are the core radius the crust thickness A/? and the total radius R as functions of L 
for a fixed total mass of O.4M0, l.OM© and \ AMiy„ respectively. It is seen that the /?, increases almost linearly with 
increasing L. The R, also increases with the increasing mass at a fixed L. This is because the stiffer the symmetry 
energy is, the larger the contribution of the isospin asymmetric part of the pressure will be, which makes the R, 
larger Moreover, the ISR decreases with the increasing L especially for light neutron stars, as the transition density 
decreases with the increasing L. As the thickness of the crust A/? and the core radius R, depend oppositely on L, the 
total radius R = Rt+ AR of the neutron star may show a complicated dependence on L. We stress here that this is the 
result of a competition between the repulsive nuclear pressure dominated by the symmetry energy contribution and 
the gravitational binding. Interestingly, it is often mentioned that the crust of neutron stars bears some analogy with 
the neutron-skin of heavy nuclei. However, they show completely opposite dependences on the L. Namely, the size of 
neutron-skin usually increases with the increasing L as a result of the competition between the nuclear surface tension 
and the pressure difference of neutrons and protons flT*], while the thickness of neutron star crusts decreases with the 
increasing L as a result of the competition between the nuclear pressure and the gravitational binding. 



The fractional mass and momenta of inertia of neutron star crust 

What is the crustal fraction AM /M of the total mass and how does it depend on the symmetry energy? Since the 
mass is simply the integration of the energy density, one expects the AM /M and AR/R have very similar dependences 
on L 1 36]. Shown in the right window of Fig. |4] (right panel) is the AM /M. The fractional mass of the crust decreases 
with the increasing L at a fixed total mass, and it decreases with the increasing total mass at a fixed value of L. The 



moment of inertia is determined by the distribution of the energy density. From the middle window, it is seen that the 
total moment of inertia increases with the increasing mass at a fixed value of L and increases with the increasing L 
at a fixed total mass. The dependence on L is relatively weak especially for light neutron stars. However, the crust 
contribution of the moment of inertia varies much more quickly with L. It decreases with the increasing neutron star 
mass at a fixed value of L and decreases with the increasing L at a fixed total mass. These are all consistent with the 
behaviors of the fractional mass and size of the crust. 

The crustal fraction of the moment of inertia A/// is particularly interesting as it can be inferred from observations 
of pulsar glitches, the occasional disruptions of the otherwise extremely regular pulsations from magnetized, rotating 
neutron stars. It can be expressed approximately as 



A/ _ 28;rP,/?3 (1 - 1.67.^-0.6^2) 
T ~ 3Mc2 1 



1 



-1 



(2) 



where m is the mass of baryons and ^ = GM/Rc^. This analytical formula has been verified by direct numerical 
calculations using both the full EOS and its PA |36]. Furthermore, it is indicated that there exists big differences for 
AI /I by comparing calculations using the full EOS and its PA |36], again due to the corresponding differences in the 
transition density. As it was stressed in ref. |6], the AI/I depends sensitively on the symmetry energy at sub-saturation 
densities through the P, and p,, but there is no explicit dependence upon the EOS at higher-densities. Experimentally, 
the crustal fraction of the moment of inertia has been constrained as AI/I > 0.014 from studying the glitches of the 
Vela pulsar [31]. Combining the observational constraint of AI/I > 0.014 with the upper bounds of p, = 0.065 fm^^ 
and Pf — 0.26 MeV/fm^ inferred from heavy-ion reactions, we can obtain a minimum radius of R > A.l + A.QM/Mq 
km for the Vela pulsar. According to this constraint, the radius of the Vela pulsar is predicted to exceed 10.5 km 
should it have a mass of 1 AMq. We notice that a constraint of /? > 3.6 + 3.9M/M0 km for this pulsar has previously 
been derived in Ref. |31] by using p, = 0.075 fm^^ and P, = 0.65 MeV/fm^. However, the constraint obtained here 
using data from both the terrestrial laboratory experiments and astrophysical observations is significantly different and 
actually it is more stringent. 



IMPRINTS OF SYMMETRY ENERGY ON GRAVITATION WAVES 

Gravitational waves are tiny disturbances in space-time and are a fundamental, although not yet directly confirmed, pre- 
diction of General Relativity. They can be triggered in cataclysmic events involving (compact) stars and/or black holes. 
They could even have been produced during the very early Universe, well before any stars had been formed, merely as 
a consequence of the dynamics and expansion of the Universe. Because gravity interacts extremely weakly with mat- 
ter, gravitational waves would carry a genuine picture of their sources and thus provide undisturbed information that 
no other messenger can deliver |50]. Gravitational wave astrophysics would open an entirely new non-electromagnetic 
window making it possible to probe physics that is hidden or dark to current electromagnetic observations |51]. 

Deformed pulsars and various oscillation modes of spherical neutron stars are among the possible sources of 
gravitational waves. In particular, the deformed pulsars are major candidates for sources of continuous gravitational 
waves in the frequency bandwidth of the existing ground-based laser interferometric detectors including the LIGO ll52ll 
and VIRGO (e.g. Ref. [53]). While the oscillations of neutron stars are mostly at frequencies much higher than 
the currently existing and planned gravitational wave observatories, their studies are of fundamental theoretical 
interest. The strain-amplitude of gravitational waves from deformed pulsars depends on the star's quadrupole moment 
determined by the EOS of neutron-rich nuclear matter On the other hand, both the frequency and the decay rate of 
the fundamental oscillation modes are determined also by the EOS of neutron-rich nuclear matter. In the following we 
present several ex amp les illustrating the imprints of the symmetry energy on gravitation waves. More details can by 
found inrefs. SmtlM]. 



Gravitational waves from deformed pulsars 

The strain amplitude of gravitational waves at the Earth's vicinity (assuming an optimal orientation of the rotation 
axis with respect to the observer) from deformed pulsars can be written as Ii58i1 

ho^ — i , (3) 
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FIGURE 5. Neutron star moment of inertia (left), quadrupole moment (middle) and ellipticity (right) as functions of the neutron 
star mass. Taken from Ref. jsiLISTll 



where v is the neutron star rotational frequency, its principal moment of inertia, e = [Ixx — Iyy)/Izz its equatorial 
ellipticity, and r its distance to Earth. The ellipticity is related to the neutron star maximum quadrupole moment (with 
m = 2) via 159D 

where for slowly rotating (and static) neutron stars <I>22 can be written as lf59|] 

In the above expression a is the breaking strain of the neutron star crust which is rather uncertain at present time and 
hes in the range a — [10^^ — 10^^] llssll . From Eqs. (O and (|4|i it is clear that Iiq does not depend on the moment of 
inertia I^^, and that the total dependence upon the EOS is carried by the quadrupole moment <I>22. Thus Eq. (O can be 
rewritten as 

, '^'22v2 

ha^X , (6) 

r 

with X — \/20457r^/15G/c^. For slowly rotating neutron stars Lattimer and Schutz ll60ll derived the following 
empirical relation for the moment of inertia 



(0.237 ± 0.008 )M/?2 



Mkm f Mkm \ 

-4.2 + 90 

MqR \MqR J 



(7) 



Using Eq. ^ to calculate the neutron star moment of inertia and Eq. ([Sj the coiTesponding quadrupole moment, the 
ellipticity e can be readily computed (via Eq. dU). The results are shown in Fig. |5] It is clearly seen that the EOS, 
especially the symmetry energy, has a strong effects on all of these important quantities. 

Since the global properties of spinning neutron stars (in particular the moment of inertia) remain approximately 
constant for rotating configurations at frequencies up to ^ 300//z |46], the above formalism can be readily employed 
to estimate the gravitational wave strain amplitude, provided one knows the exact rotational frequency and distance to 
Earth, and that the frequency is relatively low (below ^ 300//z). 

Shown in the left window of Fig. |6]is the GW strain amplitude, ho, as a function of frequency for several slowly 
rotating near-earth neutron stars. At higher rotational frequencies, relativistic rotational models have to be used ifssll . 
The solid line represents the designed upper detection limit of LIGO. Shown in the right window is the ho for the PSR 
J1748-2446 rotating at 716z assuming it has an ellipticity of £ = 10^^. In both cases clear imprints of the symmetry 
energy are seen. 



Gravitational waves from the axial w-mode oscillations of neutron stars 



In the framework of general relativity, gravitational radiation damps out the neutron star oscillations which leads 
to the frequency of the non-radial oscillations to become "quasi-normal" (complex) with a real part representing 
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FIGURE 6. Gravitational-wave strain amplitude as a function of frequency for several slow pulsars (left) and PSR J 1748-2446 
(right). Taken from Ref. L54 i55l1 




FIGURE 7. Frequency (left) and decay time (right) of w/-mode (upper) and w//-modes (lower) as functions of the neutron star 
mass M. Taken from ref. ISQI 



the actual frequency of the oscillation and an imaginary part representing the losses due to its damping 1.6 1 1 . The 
eigen-frequencies of the quasi-normal modes could be found by solving the equations which describe the non-radial 
perturbations of a static neutron star in general relativity. The critical input to solve the equation is the nuclear EOS. The 
so called w-mode associated with the space-time perturbation only exists in general relativity. It is very important for 
astrophysical applications since it is related to the space-time curvature and exists for all relativistic stars, including 
black holes. Using the MDI EOS, Wen et al. recently studied the axial w-modes fsS]. Figure [T] displays both the 
frequency and damping time of the w/- (upper frame) and w//-modes (lower frame) respectively, as a function of the 
neutron star mass. The results establish the relationship between the expected frequencies of the axial w-modes, for 
a given EOS, and the stellar mass. It is interesting to notice that there is a clear difference between the frequencies 
calculated with the MDI {x = 0) EOS and those with the MDI (x = — 1) EOS. Since the major difference between 
these two cases is the density dependence of the nuclear symmetry energy, it is obvious that the symmetry energy has 
a clear imprint on the frequencies. 

Shown in Fig. |8]are the real (upper panel) and imaginary (lower panel) parts of the eigen-frequency of wi and wjj- 
modes scaled by the mass M as a function of the compactness parameter M//?, respectively. These results suggest that 
the scaled eigen-frequency exhibits a universal behavior as a function of the compactness parameter independent of 
the EOS used. As discussed by Andersson [62], Benhar ll63ll and Tsui ll64ll this finding could be used to constrain 
the frequency and damping time of gravitational waves. This is very important for guiding the gravitational wave 
search provided the mass and radius of the prospective source (neutron star) are known. On the other hand, when 
the gravitational wave astronomy becomes a reality, namely if both the frequency and dumping time for a given 
neutron star are known this could provide information on the neutron star mass and radius. In Figs.|8]the fitting curves 
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FIGURE 8. The scaled frequency (left) and decay time (right) of vv/-mode (upper) and w//-modes (lower) as functions of the 
neutron star compactness M/R. Taken from ref. 

representing the results of Tsui et al. are also shown. It is seen that the numerical results of Wen et al. are in good 
agreement with those of Tsui et al. Ii64il . 



SUMMARY 

In summary, important progress has been made in recent years in constraining the symmetry energy with heavy-ion 
collisions. Their implications in some astrophysical phenomena have been explored. Nevertheless, the field is still at 
its beginning. While a number of potentially useful probes of the Esym{p) have been proposed, available experimental 
data are mostly for reactions with stable beams |5]. Coming experiments with more neutron-rich nuclei at several 
advanced radioactive beam facilities are expected to improve the situation dramatically. Thus, more exciting times are 
yet to come. 
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